
OKHEP-05-01 

The University of Oklahoma hep-ph/0601004 

December 2005 



Detecting Higgs Bosons with Muons at Hadron Colliders 

Chung Kao* and Yili Wangt 

Homer L. Dodge Department of Physics and Astronomy 
University of Oklahoma 
Norman, Oklahoma 73019, USA 



Abstract 

We investigate the prospects for the discovery of neutral Higgs bosons with a pair of muons by 
direct searches at the CERN Large Hadron Collider (LHC) as well as by indirect searches in the 
rare decay B s — > n + u~ at the Fermilab Tevatron and the LHC. Promising results are found for 
the minimal supersymmetric standard model, the minimal supergravity (mSUGRA) model, and 
supergravity models with non-universal Higgs masses (NUHM SUGRA). For tan/3 ~ 50, we find 
that (i) the contours for a branching fraction of B{B S — > u + u ) = 1 x 10 -8 in the parameter space 
are very close to the 5<r contours for pp — ► be/) — > 6// + /i~ + X, 0° = hP, H°, A at the LHC with an 
integrated luminosity (L) of 30 fb _1 , (ii) the regions covered by B(B S — > > 5 x 10 -9 and the 

discovery region for b(p° — > with 300 fb _1 are complementary in the mSUGRA parameter 

space, (iii) in NUHM SUGRA models, a discovery of B(B S -> ~ 5 x 10~ 9 at the LHC will 

cover regions of the parameter space beyond the direct search for befp — > bfi + fi~ with L = 300 fb _1 . 
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t E-mail address: Yili@physics.ou.edu (Y. Wang) 
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I. INTRODUCTION 



The minimal supersymmetric standard model (MSSM) [1] has two Higgs doublets <pi 
and <p2 that couple to fermions with weak isospin — 1/2 and +1/2 respectively |2|. After 
spontaneous symmetry breaking, there remain five physical Higgs bosons: a pair of singly 
charged Higgs bosons H^, two neutral CP-even scalars H° (heavier) and h° (lighter), and 
a neutral CP-odd pseudoscalar A . The Higgs potential is constrained by supersymmetry 
(SUSY) such that all tree-level Higgs boson masses and couplings are determined by just 
two independent parameters, commonly chosen to be the mass of the CP-odd pseudoscalar 
(mi) and the ratio of vacuum expectation values of neutral Higgs fields (tan/3 = v 2 /vi). 

In the MSSM and two Higgs doublet models with Model II of Yukawa interactions, the 
couplings of down type quarks and leptons with the neutral Higgs bosons are proportional 
to 1/ cos (3. Thus a large value of tan (3 greatly enhances the production rate of Higgs bosons 
produced in association with bottom quarks as well as the branching fraction of the rare 
decay B s — > fi + p~ mediated by neutral Higgs bosons. 

If tan/3 > 10, the MSSM neutral Higgs bosons are dominantly produced from bottom 
quark fusion bb -> cp° @, fl H at the CERN Large Hadron Collider (LHC). For a 
Higgs boson produced along with one bottom quark at high transverse momentum (pt), the 
leading-order subprocess is bg — > b(p° (a 0, [l(J, El, HI- If two high pt bottom quarks are 



re quir ed in association with a Higgs boson, the leading order subprocess should be gg 
pi. Il3l. E3 . Eil . 3]. We note that the importance of the process with a bottom quark was 
suggested by the authors of Ref. 13 ■ 

The LHC has a great potential to discover the inclusive muon pair channel for neutral 
Higgs bosons of minimal supersymmetry 3; 3- Recently, it was found that the discov- 



ery channel with one energetic bottom quark |2(| greatly improves the discovery potential of 
the LHC beyond the inclusive channel without bottom quarks [TtJ (pp — > 4>° — > + X) 

and the associated channel with two bottom quarks 13|| (pp — >• 660° — >• bbji + n~ + X). Since 



the ATLAS and the CMS both have excellent muon mass resolution, this discovery channel 
will provide a good opportunity to reconstruct the Higgs boson mass at the LHC with high 
precision. 

We follow the strategies developed in Ref. [2^ to investigate the discovery at the LHC of 
a neutral Higgs boson produced with one bottom quark followed by Higgs decay into a muon 
pair. We work within the framework of the minimal supersymmetric model, the minimal 
supergravity unified model, and supergravity unified models with non-universal Higgs boson 
masses at the grand unified scale (Mqut)- 



In the minimal supergravity unified model 2l( , the significance of pp —>■ <ft° —>■ + X 

is greatly improved by a large tan (3 0| because the large bb(f>° couplings make rriA and 
ran small through the evolution of renormalization group equations 23] . Consequently, the 
production cross section is further enhanced by a large value of tan (3. 
The rare decay B s — > has a small branching fraction 

B(B S -> p + fx~) = 3.4 x 1(T 9 (1) 

in the Standard Model (SM) of electroweak interactions 0,|2^. A recent calculation [26| 

suggests 

B(B S -> n + fiT) = (5.1 ± 1.1) x HT 9 (2) 
with updated parameters. The current experimental upper limit is 

B(B S -> p + p-) < 1.5 x 1(T 7 (3) 
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obtained by the Collider Detector at Fermilab (CDF) and the D0 collaborations 27(. While 
this branching fraction is small in the SM, it could become large in supersymmetric models 
28] -[46] and this rare decay provides a possible opportunity for the CDF and the DO 



experiments to discover new physics in the near future. 

In this letter, we investigate the discovery potential of the direct searches for the Higgs 
bosons pp — > bcjP — > bfi + fi~ +X at the LHC and that of the indirect searches for Higgs bosons 
in B s — > fi + fi~ at the Fermilab Tevatron Run II within the framework of supersymmetric 
models. We make three major contributions: (a) studying the LHC discovery potential for 
pp — > b<f)° — > b[i + [i~ + X in the minimal supergravity model (mSUGRA) and in supergravity 
models with non-universal Higgs bosons masses at Mgut (NUHM SUGRA), (b) evaluating 
B(B S — > in NUHM SUGRA models, and (c) comparing these two promising channels 

in the MSSM, the mSUGRA, and non-universal SUGRA models. In Section II, we discuss 
our strategies and results for the minimal supersymmetric standard model. Sections III 
presents the discovery contours in the parameter space of the mSUGRA model as well as 
that of the NUHM SUGRA models. Conclusions are drawn in Section IV. 



II. THE MINIMAL SUPERSYMMETRIC STANDARD MODEL 

In this section, we consider the direct searches for pp — > b<f)° — > b[i + fi~ + X at the LHC 
and indirect searches for Higgs bosons in B s — > at the Fermilab Tevatron Run II 

within the framework of the minimal supersymmetric standard model. 



A. b<p° -> bfi+fj,- 

Applying previous calculations [20] for the Higgs signal at the LHC we evaluate the cross 
section of pp — > b(f)° — > bfi + fi~~ +X with the Higgs production cross section a{pp — > b(f)° + X) 
multiplied by the branching fraction of the Higgs decay into muon pairs B(<f)° — > 
The cross section for pp — > 60° + X (0° = H °, h°, A ) via bg — > 60° is calculated with 
the parton distribution functions of CTEQ6L1 |4jj with a factorization scale [if = Mh/4: 
0, Hij]. Unless explicitly specified, b represent a bottom quark (b) or an anti-bottom quark 
(b). We have also taken the renormalization scale to be M#/4 that effectively reproduces 
the effects of next-to-leading order (NLO) ^(j with a K factor of one for the Higgs signal. 
The bottom quark mass in the cj)°bb Yukawa coupling is chosen to be the NLO running mass 
m b (fifi) [49], which is calculated with m b (pole) = 4.7 GeV and the NLO evolution of the 
strong coupling |50j. 

In our analysis, we consider dominant physics backgrounds to the final state of bfi + fi~ 
from bg — > b/i + }j,~ (bfifi) as well as gg — > bbW + W~ and qq — > bbW + W~ (bbWW) followed 
by the decays of W ± — > y^v^. In addition, we have included the background from bg — > 
bix^vyrv and bg — > bfi~vfi + b>, which has major contributions from bg —>■ tW~ and bg — > iW + 
{tW). The muons from b decays can be removed effectively with isolation cuts [13]. We apply 
a K factor of 1.3 for the b\x\x background j5l|, a K factor of 2 for bbWW 5^, 53|, and a K factor 



of 1.5 for tW |54j. Furthermore, we consider backgrounds from pp — > jfi +X,j = g,q or 
q with q = u,d, s, c, where a jet is mistagged as a b with a K factor of 1.3 for these processes. 

We adopt the acce pta nce cuts as well as the 6-tagging and mistagging efficiencies of the 
ATLAS collaboration [19(. In each event, two isolated muons are required to have Pt(^) > 20 
GeV and 1 77 (/x) | < 2.5. 
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For an integrated luminosity (L) of 30 fb -1 , we require (i) pr(b,j) > 15 GeV, (ii) 
\r)(b,j)\ < 2.5, and (iii) the missing transverse energy ( lfi T ) should be less than 20 GeV 
to reduce the background from from bbWW and tW which contains neutrinos. The b- 
tagging efficiency (e?,) is taken to be 60%; the probability that a c-jet is mistagged as a 6-jet 
(e c ) is 10% and the probability that any other jet is mistagged as a 6-jet (e,,) is taken to be 
1%. For < 100 GeV, we change the requirement to pr(fJ') > 10 GeV for muons in both 
the Higgs signal and the background. 

For a higher integrated luminosity of 300 fb -1 , we require pr(b, j) > 30 GeV and e& = 50%. 
In addition, the missing transverse energy ( ]ft T ) in each event should be less than 40 GeV. 

To study the discovery potential of pp —>■ b<f)° + X —>■ b[i + fi~ + X at the LHC, we 
calculate the background from the SM processes of pp — > bp + u~ + X in the mass window 
of ± A¥ M+r where AM A+r = 1.64[(r^/2.36) 2 + o^] 1 / 2 0. T is the total width of 
the Higgs boson, and a m is the muon mass resolution which we take to be 2% of the Higgs 
boson mass |19I ]. 



B. B s — > n+fj,- 

In our analysis for B s — > U + U~ within the framework of minimal supersymmetry, we 
follow the approach in Refs. ji^, 42|| and adopt the formulas in Ref. [42]. We make the 



following choices: (i) The matrix of the up-type Yukawa couplings is diagonal, (ii) The 
down-type Yukawa coupling matrix is Fp = DVq KM , where D is the matrix diagonalized 
from Fb and Vckm is the Cabibbo-Kobayashi-Maskawa matrix, (iii) We neglect the masses 
of the d and the s quarks as well as terms that are second order or higher in the Wolfenstein 
parameter A. (iv) At the tree level, the CKM matrix is the only source for flavor changing 
neutral current (FCNC), and (v) We include FCNC contributions from one-loop diagrams 
involving charginos and up-type squarks as well as gluino and down-type squarks. 

In addition, we adopt a common mass scale for supersymmetric (SUSY) particles and 
parameters Msusy = Trig — — p, — —Af, where Af = A t = = A T are the trilinear 
couplings for the third generation. Two values of Msusy are considered: (a) Msusy = 350 
GeV and (b) M SUS y = 1000 GeV. 

In Figure 1, we present the contours for the branching fraction in the MSSM B(B S — > 
/i + /i _ ) = 1.5 x 10~ 7 (current experimental limit), 3 x 10~ 8 , 1 x 10~ 8 , and 5 x 10~ 9 as well as 
the discovery contours of b(p° — > bji + fi~ for integrated luminosities of 30 fb -1 and 300 fb -1 
at the LHC in the (m^, tan (3) plane for two values of common masses: (a) Msusy — m g — 
m f - = 350 GeV, and (b) M SUS y = m- g = m f ~ = 1000 GeV. 

We note that for M S usy = 350 GeV, the LHC will be able to discover pp — > b(f)° — > 
bfi + {i~ + X with an integrated luminosity (L) of 30 fb -1 in a significantly large region of 
the parameter plane beyond B(B S — > = 3 x 10~ 8 . If the gluino and scalar fermions 

have a common mass of approximately 1 TeV then the contour for a branching fraction 
of B(B S — > /U + /i _ ) = 1 x 10 -8 in the parameter plane is very close to the 5a contour for 
pp -> b<f b{i + n~ + X at the LHC with L = 30 fb -1 . 

Furthermore, with a higher luminosity of 300 fb -1 , the LHC will be able to discover 
pp — > b<f)° — > bfi + fi~ + X for Msusy = 1000 GeV in a very large region of the (m^, tan ft) 
plane. The discover contour for high luminosity with a large Msusy is very close to the 
contour for B(B S — > /U + /i - ) = 5 x 10~ 9 that is the not far away from the SM expectation. 
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FIG. 1: Discovery contours for pp — > b(jr — > b/j,p,+ X at the LHC and contours of the branching of 
-B s — > /U + /*t~ in the minimal super symmetric standard model for (a) m§ = mj = 350 GeV = —^4/ 
and (b) nig = mj = 1000 GeV = —Af. The discovery region is the part of the parameter space 
above the contour. 



III. SUPERSYMMETRIC UNIFIED MODELS 



In this section, we consider both pp — > b(f)° — > + X and B s — ► in the minimal 

supergravity model and supergravity models with non-universal boundary conditions for the 
Higgs boson masses at the grand unified scale (Mqut)- We evolve supersymmetry masses 
and couplings from the grand unified scale using two-loop renormalization equations in 
ISAJET 7.72 |55[ to calculate MSSM masses, mixing angles and couplings. The following 
theoretical requirements are imposed on the evolution of renormalization group equations: 
(i) radiative electroweak symmetry breaking (EWSB) is achieved; (ii) the correct vacuum for 
EWSB is obtained (tachyon free); and (iii) the lightest SUSY particle (LSP) is the lightest 
neutralino (x?)- 



A. The minimal supergravity unified model 

In the minimal supergravity (mSUGRA) model j2jj , supersymmetry is broken in a hidden 
sector and SUSY breaking is communicated to the observable sector through gravitational 
interactions. The mSUGRA parameters are chosen to be a scalar mass (mo), a gaugino mass 
( m i/2), a trilinear coupling (A ), the sign of a Higgs mixing parameter (//), and the ratio of 
Higgs field vacuum expectation values at the electroweak scale (tan/3 = Vijv\). The value 
of A Q only significantly affects results for high tan/3; we initially take A = and study the 
A dependence later. 

Figure 2 displays the discovery contours of b<p — ► bfi + [i~ for an integrated luminosity of 
30 fb~ x and 300 fb _1 at the LHC as well as contours for four values of the branching fraction 
B(B S — > /U + /i~) in the (m^, mo) plane of the mSUGRA model for four values of tan f3 = 20, 
30, 40, and 50. Also shown are the parts of the parameter space (i) disfavored by theoretical 
requirements or (ii) excluded by the chargino search at LEP 2 with m x ± < 103 GeV. 

There are several interesting aspects to note in Figure 2. 

(i) If tan/? < 30, only a tiny region of the parameter space with small values of m 1( / 2 and 
mo will likely lead to observable signals for either B s — ► at the Tevatron Run II 



5 



mSUGRA 




FIG. 2: Discovery contours for pp — > bcjr — > + X at the LHC and contours of the branching 
of B s — > in the minimal supergravity unified model for (a) tan/? = 20, (b) tan/3 = 30, (c) 

tan/3 = 40, and (d) tan/3 = 50, Also shown are the parts of the parameter space (i) excluded by 
theoretical requirements (slant-hatched and dark shaded), or (ii) excluded by the chargino search 
at LEP 2 (horizontally-hatched). 



or b(f)° — » bjj, + fi at the LHC. 

(ii) For tan f3 < 40, direct searches for bcfP — ► bfi + [i~ at the LHC with L = 30 fb _1 covers a 
much large region in the mSUGRA parameter space than B(B S — > /U + /i~) > 1 x 10~ 8 . 

(iii) If tan/3 50, the discovery contour for b<p° — > b[i + [i~ at the LHC with L = 30 fb _1 is 
very close to the contour for B(B S — > fi + n~) > 1 x 10 -8 in the mSUGRA parameter 
space. In addition, both discovery channels at the LHC become complementary. The 
direct searches for b(f>° — > b[i + fi~ with L = 300 fb _1 covers a significant region beyond 
the contour of B(B S — > /i + /i~) = 5 x 10 -9 . Likewise, the rare decay with B(B S — ► 
H + fj>~) > 5 x 10~ 9 covers a large region beyond the discovery contour of the direct 
search for b<ft° — > bfi + {i~ with L = 300 fb _1 . 



B. The mSUGRA model with non-universal Higgs masses 



In our analysis for non-universal supergravity models, the GUT-scale Higgs masses are 
parameterized as 

m% l (GVT) = (l + 6 i )m%, z = l,2. (4) 

The nonuniversality of Higgs-boson masses at Mqtjt can significantly affect the values of 
Higgs masses and couplings at the weak scale 5(| 53, HI 59l |. 

We find that a decrease in with a negative Si as well as an increase in m# 2 with a 
positive S2 at Mqut will lead to a smaller mass at the electroweak scale for the Higgs pseu- 
doscalar (A ) or the heavier Higgs scalar (H°) than that in the mSUGRA model. There- 
fore, we choose three sets of values for Si to study the discovery potential for detecting 
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Higgs bosons with muons in SUGRA models with non-universal Higgs boson masses: (i) 
S 1 = -0.5, 8 2 = 0, (ii) 5i = 0, 8 2 = 0.5, and (hi) Si = -0.5, 5 2 = 0.5. 

In Figure 3, we present the discovery contours of b<f)° — > bfi + fi~ for integrated luminosities 
of 30 fb -1 and 300 fb -1 at the LHC as well as contours for four values of the branching fraction 
B(B S — > in the (777.1/2, m o) plane for a NUHM SUGRA model with S\ = —0.5 and 

62 = with tan/3 = 20, 30, 40, and 50. In addition, we show the regions of the parameter 
space (i) disfavored by theoretical requirements or (ii) excluded by the chargino search at 
LEP 2 with m + < 103 GeV. 




m„ 2 (GeV) m 1/2 (GeV) 



FIG. 3: The 5a contours for pp — » b(jr — > bfifi + X at the LHC with an integrated luminosity of 30 
fb -1 and 300 fb -1 as well as contours for the branching fraction of B s — » in the (mirt, ttiq) 

plane of a non-universal SUGRA model with fi > 0, A$ = and non-universal boundary conditions 
5 1 = -0.5 and 5 2 = 0, for (a) tan (3 = 20, (b) tan/3 = 30, (c) tan/3 = 40, and (d) tan/? = 50, Also 
shown are the parts of the parameter space (i) excluded by theoretical requirements (slant-hatched 
and dark shaded), or (ii) excluded by the chargino search at LEP 2 (horizontally-hatched). 

Figure 4 shows contours for four values of the branching fraction B(B S — > n + fi~) in the 
(777.1/2,7770) plane in a supergravity unified model with S± — and 62 = 0.5. Also shown are 
the discovery contours of pp —>■ b<f)° — * bfi + fi~ + X for integrated luminosities of 30 fb -1 and 
300 fb -1 at the LHC in the (7711/2, mo) plane for four values of tan/3 = 20, 30, 40, and 50. 

If both Higgs boson masses are different from the common scalar mass at Mgut ; then 
theoretically favored region shrinks greatly. We present contours for four values of the 
branching fraction B(B S — > /i + /U~) in the (7711/2,7710) plane in a supergravity unified model 
with 5\ = —0.5 and 62 = 0.5 in Figure 5. In addition, we show the discovery contours of 
pp — > b(p — > bfi + fi~ + X for integrated luminosities of 30 fb -1 and 300 fb -1 at the LHC in 
the (777.1/2,777,0) plane for four values of tan/3 = 20, 30, 40, and 50. 

In all three NUHM SUGRA cases that we have considered, itia and m# are smaller 
than those in the mSUGRA model for the same values of mo and 777.1/2. Consequently, both 
b<f)° — > b{i + fi~ and B s — > fi + fi~ will be able to cover regions of the parameter space with larger 
values of m and 777,1/2. We note that for tan/3 > 50, the observable region for b(f>° — > 
at the LHC with L = 30 fb -1 is comparable to that of B(B S — > /U + /i~) > 1 x 10~ 8 . However, 
the contour for B(B S — > /i + /i _ ) = 5 x 10 -9 lies beyond the discovery contour for b(f>° — > /i~ 
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FIG. 4: The ha contours for pp — > bcfP — > bfifl + X at the LHC with an integrated luminosity of 30 
fb _1 and 300 fb -1 as well as contours for the branching fraction of B s — > in the (m^ 2 ,mo) 

plane of a non-universal SUGRA model with \x > 0, Aq = and non-universal boundary conditions 
St = and 5 2 = 0.5, for (a) tan/3 = 20, (b) tan/3 = 30, (c) tan/3 = 40, and (d) tan/? = 50, Also 
shown are the parts of the parameter space (i) excluded by theoretical requirements (slant-hatched 
and dark shaded), or (ii) excluded by the chargino search at LEP 2 (horizontally-hatched). 



at the LHC with L = 300 fb _1 . 

In the NUHM SUGRA model with <5i = -0.5 and S 2 = 0, (m Hl = 0.707m and m H2 = 
mo), most of the (mi/ 2 ,m ) plane is theoretically favored for tan/3 < 40. If m# 2 is larger 
than m with <5 2 = 0.5, the theoretically disfavored region grows rapidly as the value of tan (3 
increases. 



IV. CONCLUSIONS 

In supersymmetric models, the muon pair discovery channels offer great promise for the 
detection of indirect Higgs signatures in B s — > [i + at the Fermilab Tevatron as well as for 
the direct signal of pp — > 60° — > + X at the CERN LHC. 

If scalar fermions and gluino have a mass close to 1000 GeV, then the exclusion contours of 
the Tevatron search for B(B S — > ^ + ^) — 1 x 10~ 8 are comparable to the discovery contours 
of the LHC search for 60° — > b[i + pT with an integrated luminosity of 30 fb -1 . However, if 
SUSY particles have a mass close to 350 GeV, the direct search for 60° — > bp + p~ at the 
LHC becomes much more promising than B s — > p + p~ . 

In supergravity unified models, the branching fraction of B s — > p + p~ and the significance 
of pp — > b<ft° —>■ bfi + fi~ + X are greatly improved by a large tan j3 because the large bb<p 
couplings make tjia and m# small through the evolution of renormalization group equations 
and enhance the production cross section for Higgs bosons even more. In the mSUGRA 
model and in supergravity models with non-universal Higgs masses, the direct signal of 
b<f)° — > bp^ pT at the LHC can be discovered with a luminosity of 30 fb -1 in a large region of 
the parameter space that is comparable to that of B{B S — > p + p~) = 1 x 10~ 8 for tan (3 < 50. 

For a large value of tan /?, the Tevatron Run II with an integrated luminosity of 10 fb -1 
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FIG. 5: The 5a contours for pp — » 60° — > + X at the LHC with an integrated luminosity of 30 
fb _1 and 300 fb -1 as well as contours for the branching fraction of B s — » in the (m^ 2 ,mo) 

plane of a non-universal SUGRA model with // > 0, Aq = and non-universal boundary conditions 
6 1 = -0.5 and 5 2 = 0.5, for (a) tan (3 = 20, (b) tan (3 = 30, (c) tan (3 = 40, and (d) tan (3 = 50, Also 
shown are the parts of the parameter space (i) excluded by theoretical requirements (slant-hatched 
and dark shaded), or (ii) excluded by the chargino search at LEP 2 (horizontally-hatched). 



will be able to observe the indirect Higgs signal with B(B S — > /U + /i~) ^ 1 x 10~ 7 j45( which 
then can be confirmed by the direct search of pp — > b<p° — > b^i + + X at the LHC [20] . 
Even if there are no signs of new physics at the Tevatron Run II, we will be able to set 
meaningful limits on important parameters such as the Higgs masses and the ratio of the 



vacuum expectation values of Higgs fields tan/? = v%lv\ |43| for the MSSM with minimal 
flavor violation. The Tevatron Run II with an integrated luminosity of 8 fb _1 will be able 
to exclude B[B S — > ^ 3 x 10~ 8 [6(} which then can provide important guidance to 

detect pp -> 60° -> 6^+^ + x at the LHC 0. 

If tan /3 ^ 50, the regions covered by B(B S — > /i + /i _ ) > 5 x 10~ 9 and the discovery region 
for 60° — > b{i + with 300 fb _1 are complementary in the mSUGRA parameter space: the 
direct searches for h(f)° — > bfi + fi~ with L = 300 fb _1 can cover a significantly large region 
beyond B(B S — > /i + /i _ ) > 5 x 10~ 9 , and vice versa. However, in NUHM SUGRA models, 
a discovery of B(B S — > fi + n~) — 5 x 10~ 9 at the LHC will cover regions of the parameter 
space beyond the direct search for 60° — > bfi + fi~ with L = 300 fb _1 . 
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